Cosmic Microwave Background Delensing Revisited: Residual Biases and a Simple Fix 
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The delensing procedure is an effective tool for removing lensing-induced B-mode polarization 
in the Cosmic Microwave Background to allow for deep searches of primordial B-modes. However, 
the delensing algorithm existing in the literature breaks down if the target B-mode signals overlap 
significantly with the lensing B-mode (~ 300 < I < 2000) in multipole-Z space. In this paper, we 
identify the cause of the breakdown to be correlations between the input B-map and the deflection 
field estimator (EB). The amplitude of this bias is quantified by numerical simulations and com- 
pared to the analytically derived functional form. We also propose a revised delensing procedure 
that circumvents the bias. While the newly identified bias does not affect the search of degree scale 
B-mode generated by tensor perturbations, the modified delensing algorithm makes it possible to 
perform deep searches of high-i B-modes such as those generated by patchy reionization, cosmic 
strings, or rotation of polarization angles. Finally, we estimate how well future polarization experi- 
ments can do in detecting tensor- and cosmic string- generated B-mode after delensing and comment 
on different survey strategies. 



I. INTRODUCTION 

Linear, scalar perturbations do not generate B-mode 
polarization in the cosmic microwave background (CMB) 
[l| @ • This exact symmetry makes B-mode polarization 
a power tool to look for ingredients beyond standard 
cosmology, such as tensor perturbations (gravitational 
waves) , cosmic strings, patchy reionization, and polar- 
ization rotation Q [3| @ 0- Gravitational lensing 
by large scale structures converts B-mode to B-mode, 
acting as a foreground contamination. Because the fre- 
quency dependence of lensing B-mode is that of a 2.73K 
blackbody, it cannot be removed with multi-frequency 
observations and will eventually limit future B-mode ob- 
servations [|[ 

Fortunately, because of the unique non-Gaussian sig- 
natures in the lensing B-mode, it can be separated from 
the primordial Gaussian B-mode [H| ■ This "delens- 
ing" procedure reconstructs the intervening lens by look- 
ing at the higher order statistics in the temperature and 
polarization maps [HJ [3 0, calculates the expected 
lensing B-mode from the observed B-mode and the re- 
constructed mass distribution, and subtracts it from the 
observed B-mode [15J. It has been shown that this pro- 
cedure is effective in removing the lensing B-mode [la ]. 
External data sets can be used to gain information on 
the matter distribution along the line of sight, but if the 
instrument noise is low (<~ 3/iKarcmin) [rn ], internal 
delensing from the CMB becomes much more efficient 



than delensing with external data. 

In this paper, we point out that because of higher or- 
der effects that have been ignored in previous treatments, 
the delensing algorithm as described in the literature fails 
if the target B-mode signals reside in the same I range 
as the lensing B-mode (~ 300 < Z < 2000). This does 
not create a major problem for the search of degree scale 
B-mode generated by tensor perturbations, because the 
observed B-mode map can be high-pass filtered (above 
I ~ 150) before the lensing reconstruction 0]. How- 
ever, the current delensing algorithm is not applicable 
to the search of high-Z B-modes such as those generated 
by patchy reionization, cosmic strings, or rotation of po- 
larization angles. In this paper, we identify the cause of 
the breakdown to be a sizable correlation between the in- 
put B-map and the deflection field EB estimator. This 
is important because at the sensitivity level achievable 
by next generation ground based experiments, the most 
sensitive estimator for the deflection field is the EB esti- 
mator 0]. We explicitly derive the functional forms for 
a few of the bias terms and verify them with numerical 
simulations. The phenomenology caused by the bias is 
also explored and discussed. 

In section [V] we propose a revised delensing procedure 
that takes cares of all the high order terms that we have 
not analytically calculated, and show in simulations that 
it avoids the bias. This method is straightforward to 
implement, and manifestly bias-free. Furthermore, the 
new algorithm can be applied to quadratic estimators 



2 



- Lensed B power 

- String B power (G|i=10~ 7 ; 

Tensor B power (r=0.05) 
Noise (A =1.41|iK-arcmin,FWHM=4 arcminV, 




10" 



C\l 
T3 _ 

o 10 



10" 



1(1+1)0^/231 

+ Quadratic (simulation) 
- - - Quadratic (theory) 

Iterative (simulation) 
— ■ Fisher limit 

++ + +'++_++ +++ 

■ - - - ~ ; t M-±±+ " - 



500 



1000 



1500 



FIG. 1. This figure shows the shapes of B-mode power spec- 
tra generated by tensor perturbation (tensor-to-scalar ratio 
r = 0.05), cosmic strings (strings tension Gfi = 10 -7 ) and 
the gravitational lensing. While the tensor-induced B-mode 
and lensing B-mode are well separated in I space, the lensing 
B-mode overlaps with B-mode from strings and cannot be 
delensed using high-pass filters. The same is true for rota- 
tion generated B-mode [f| , which would have an B-mode like 
spectrum that extends to high-Z (not shown). 



FIG. 2. The comparison of simulated reconstruction error 
(or bias on power spectrum of deflection field) with theory 
using quadratic and iterative estimators (EE and EB). The 
theoretical estimate of error and Fisher limit are from Eq. 
(66) and (32) in [14j. The instrument noise level is set to 
be A p — 1.41/xK arcmin in these simulations, with FWHM= 
4'. Here we use 64 realizations of CMB fields but only one 
realization of lensing potential. We choose the convergence 
factor to be 0.12 and perform 64 iterations, as suggested in 




as well as likelihood-based lensing estimators (iterative 
estimators), allowing the lowest possible residuals to be 
achieved. In section I VII we apply the results obtained 
in section |V| to forecast the sensitivity of future exper- 
iments to tensor-induced and cosmic string-induced B- 
polarization with iterative delensing. 

Throughout this paper, We use a flat ACDM cos- 
mology with parameters flt,h 2 = 0.0226, il c h 2 = 0.11, 
h = 0.7, n s = 0.96, A^(fc ) = 2.46 x 10~ 9 at k = 0.002 
Mpc -1 , and r = 0.09. The tensor-to-scalar ratio is de- 
fined as initial tensor/scalar power spectrum amplitude: 
r = A^(fco)/A^(fco). The tensor and lensing i3-mode 
power spectra are produced by CAMB [18j, and the 
string B-mode power spectrum is using CMB ACT [TJ|. 
All simulations were carried out under the flat-sky ap- 
proximation. In this paper, we focus on the treatment 
of delensing bias and defer the complete treatment of de- 
lensing with EB leakage from map boundaries to a future 
paper. 



II. LENS RECONSTRUCTION 

In this section we provide a brief overview of gravi- 
tational lensing and deflection field reconstruction. For 
more details please refer to 13 1 [l4| (l7| . 

The large scale structures deflect the CMB photons 
randomly as they propagate through the universe. Since 



the deflection angles are small (~ 2 arcmin), we consider 
only the first order density perturbations and apply the 
Born approximation. We can then model the lensed fields 
as: 



X(n) =X(h + d(h)) 



(1) 



X(n) can be CMB temperature T or Stokes parameters 
Q and U. The tilde indicates the unlensed fields and d(h) 
is the deflection angle which is the gradient of integrated 
lensing potential 4>(n) (for more detail see 

The gravitational lensing breaks the statistical isotropy 
in the unlensed CMB. In fact, this is what makes the re- 
construction of the deflection field possible. The quantity 
V • (XVX) gives us the best quadratic estimator to de- 
tect the anisotropy caused by lensing [TJ[ , where X is the 
filtered field (C r " 



xx 



C NN )- 1 X. The quantity C 



NN 



the noise covariance matrix, and C xx is the covariance 
matrix of the lensed field averaged over lensing potential 
realizations so that it is diagonal in harmonic space. X 
is another filtered field C XX {C XX + C NN y 1 X 1 where 
C xx is the covariance matrix of unlensed field. 

A more complicated likelihood-based method that goes 
beyond the quadratic order has also been introduced to 
improve the lens reconstruction [ljj]. This is especially 
relevant for future CMB experiments that have low in- 
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strument noise (a few arcmin). The main distinc- 
tion between this approach and the quadratic estima- 
tor is that the likelihood method replaces the weighting 
{C xx + C NN )- 1 by (C xx + C NN )~ 1 and introduces 
mean field subtraction. As pointed out in [21], this also 
provides a way to deal with non-uniform sky coverage. 
The likelihood-based method can also be understood as 
an iterative process that recalculates the quadratic esti- 
mator repeatedly. The lensing-induced B-modes can be 
removed iteratively until the series converges to a resid- 
ual B-mode map that depends on instrument noise level. 

We developed a simulation pipeline based on the pro- 
cedures outlined above and tested it with a variety of in- 
put parameters. We followed the algorithms in [l3j and 
[l4j to simulate EE and EB quadratic and likelihood- 
based (iterative) estimators and obtained results consis- 
tent with those in their papers (see Fig. [5]). Through- 
out this paper, our simulations focus on noise level: 
A p = 1.41/iK arcmin, FWHM= 4' and sky coverage: 
10° x 10°, with grid spacing of 1 arcmin, unless otherwise 
stated. This noise level is achievable in future ground- 
based experiments such as the POLAR Array. It is worth 
mentioning that the analytic estimate for the quadratic 
estimator error [l3[ is lower than the simulated results. 
This is the result of approximating E and B modes as 
Gaussian, a known effect reported in [22| [23| (24|. The 
analytic estimate for the reconstruction uncertainty in 
iterative estimator is more difficult. Here we plot the 
Fisher limit to provide the lower bound of the lens re- 
construction error (see [l4| for details). 

III. DELENSING AND THE CAUSE OF THE 
DELENSING BIAS 

After reconstructing the lensing potential, or equiva- 
lently the deflection field, we can use this information 
to either numerically reverse lensing to obtain unlensed 
TQU maps, or to calculate the expected B-polarization 
induced by lensing and subtract it from the observed B- 
map. We tried both methods in simulations and did not 
find significant difference between the two at low I region, 
before the effects of pixelization kicks in. We choose the 
latter method (subtraction-based) for most of the simula- 
tions since it has a nice analytic representation in Fourier 
space. In practice, we Wiener-filter the observed E-mode 
(lensed and with instrumental noise) and approximate it 
as the unlensed E-mode. This approximation is justified 
because weak lensing only alters E-mode slightly. We 
then apply the deflection field estimator to "lens" this 
E-map to calculate the lensing induced B-mode. To first 
order in lensing potential </>, the expected B-mode from 
lensing is given by 

[1) J (2tt)2 1 {ll ' l2> C EE + N EE cf* + Nf* ' 

(2) 



Nf^ is lens reconstruction error. For EB estimator 
and approximating E and B modes as Gaussian we can 
express it as: 

n m = \[ th fJMHc^f 1 1 

1 [J (2tt) 2 (Cf 1 - + N EE^ C B lm . +N BB) 

(3) 

where fih^h) — sin2(pi t i 2 {l2 • I) (I — l± + 1%), cj> is the 
estimated lensing field, C x lens and C l lonB are the lensed 
EE and BB power spectra, ipi^ = fi l ~(pi 2 (<pi is angle 
between I and x-axis). This estimator is exactly Eq.(5) 
in [171 ] under flat sky approximation. The estimator for 
delensed B-mode is simply the difference between the 
observed and estimated B mode 

B dc \l) = B ohs (l)- B lcns {l). (4) 

The observed B-mode B ohs (l) includes B-mode from 
lensing B lcns (l), noise B nolse (/), and primordial B-mode 
B prl (l) if it exists. The primordial power spectrum based 
on squaring this estimator (Cf dcI = (B dol *(/)B dol (/))) 
will be biased. The sources of biases are imperfect lens 
reconstruction, instrument noise in the observed E, B 
maps, and their complicated cross power spectra. Out 
of the 10 quadratic terms in Cj B<ul , we know B lens (l), 
B nmsc (l) and B prl (l) do not correlate with one another. 
As a result, there are only 7 non-trivial terms. Two of 
them are auto-spectra, primordial B-mode power cf pT1 
and noise power (7 ; Bnoisc . We denote the total contribu- 
tion from the rest of the 5 terms as C^ Ioa . Explicitly, 

Cf r ™ = Cf 1 ™ 5 + (B lons *(Z)B lcns (Z)) -2(B lcns *(Z)B lons (0) 

- 2(B pri *(/)B lcns (/)) - 2(B noise *(/)B lcns (0). (5) 

All these terms will contribute to the lensing residuals, 
including the negative contributions from the cross cor- 
relations. Knowing the delensing biases for a given set 
of experimental parameters is important for planning of 
future CMB experiments. Furthermore, the properties 
of the residuals, including their correlations, need to be 
understood in actual data analysis so that significance of 
the target signals can be evaluated. An approximated 
residual for an EB quadratic estimator has been given 
analytically by [l7| : 




(6) 

where I2 — I + h- This estimated residual apparently 
came from the first three terms in Eq. ([5]), represent- 
ing imperfect lens reconstruction (JV. ) and instrument 
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FIG. 3. These plots demonstrate the delensing biases caused by high order correlations, (a): the blue line is tensor B -mode 
power (r = 0.05) and the red dot-dashed line is delensing residuals predicted by Eq. first given in [l7| . The green dashed 
line is theoretical prediction of 2(B pri * (l)B lens (l)} given by Eq. 0. The black crosses are simulation result using an EB 
quadratic estimator (400 realizations). This shows that the bias is far from negligible compared to the residual, (b): The blue 
line is lensing B power and the red dot-dashed line is noise power (A p = 1.41/xK arcmin, FWHM=4 arcmin). the green dashed 
line is theoretical prediction of 2(B noisc * (l)B lcns {l)) given by Eq. ([8]). The black crosses are simulation result using an EB 
quadratic estimator (400 realizations). Again, the bias is almost 60% of the total residual. 



noise (N t • ). As these two terms vanish, the resid- 
ual vanishes as well. In deriving this residual [17] ig- 
nores 2{BP"*(l)B lcns (l)) and 2{B noisc *(l)B lcns (l)}, even 
though they admited it has not been justified. They 
also approximated the calculation of (B lcns * (I) B lens (I)) 
and (B lcns *(l)B lcns (l)) by neglecting (E ohs E ohs <p(j)} c con- 
nected four-point function. 

The residual for iterative estimator can be evaluated 
in the same way using Eq. ([3]) and Eq. ([6]). The only 
difference is iterative estimator improves the ability of 
lens reconstruction by depressing C^' " 3 to C Blcs . There- 
fore, C El °" will appear simultaneously in the left-hand 
and right-hand sides of Eq. ([6]) . We can use an iterative 
algorithm to solve it. 



IV. THE MAGNITUDE OF DELENSING BIAS 

In this section, we show the neglected terms ( 
{BP"*(l)B lcns (l)) : {B Doisc *(l)B lcns (l)) and (E ohs E ohs <f>4>) c 
) contribute significantly to C El °" under parameters 
that are typical of future experiments. In the simplest 
case, we obtain <f> from EB quadratic estimator and 
the estimated lensed -B-mode is quadratic estimator of 
E ohs and cj>, {BP ri *(l)B lcns (l)) will be four-point function 
(BP"E ohs B ohs E ohs ). By applying Wick's theorem and 
ignoring EB cross correlation, we obtain 



(B pli *(l)B lcns (l)) = 

C B lma . +N BBj ( 27r )2 C EE + N EE {C+f + N+fy 

' (7) 

where I2 = I + h. Here we ignore the non-scalar i?-mode 
because it will be much smaller than scalar E-mode. In 
practice, the bias term contributes negatively to the de- 
lensed power Eq. ([5]) with an I dependence similar to the 
primordial B-mode. 

Since J5 nolse (Z) acts like B p "(l) in the cross correlation, 
the cross correlation (B noisc *(l)B lcns (l)) is in the form 
similar to Eq. ([7]) except that cf v " is replaced by N^ B , 

(B noisc *(l)B lcns (l)) = 

C B )m+N BBj ( 27r )2 c ee +n ee (cfy + Nfyy 

2 (8) 

Our simulation results (Fig[3|) show good consistency 
with the analytic calculation. Up to ~ 60% of pri- 
mordial _B-mode power spectrum will be supressed (at 
I - 60) by the cross correlation term -2(BP"* (l)B lcns (l)) 
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(the location of recombination peak of tensor .B-mode) . 
Even the instrumental noise is artificially "delensed" 
by a similar fraction because of the contribution from 
_2(B aoise *(l)B lens (l)) in Eq. ©. As a result, delensing 
simply fails to improve signal-to-noise ratio of the pri- 
mordial B-mode. 

The calculation of the additional bias term 
{E ohs E ohs <j)(t)) c is daunting, if is obtained by 
an EB estimator. It will be an 8-point function 
^obs^obs^obs^obs^pri^pri^^ Therefore, we choose to 

do simulation to see if 2(B nolsc *(/)B lcns (Z)) can account 
for all the biases in estimating CP ial , in the case of 
zero tensor .B-mode. In FigfJJ We follow the delensing 
procedure described in section IIIII to get (7 ; Bdcl for both 
cases with r = 0.05 and r = 0. We have removed E 
and B modes at I < 144 in lens reconstruction. We 
can see that even though Cf dal has been corrected by 
the simulated 2(B noisc *(l)B lcns (l)), the resulting Cf dcl 
is still lower than C^' as estimated in Eq. ©. This 
indicates that the contribution from (E ohs E ohs <fxf>) c is 
not neglegible. 

All these biases and correlations can be avoided for 
searches of tensor-induced -B-mode, because the modes 
located at low / can simply be removed prior to lens- 
ing reconstruction. As a result, ij> contains no B obs and 
B ohs at low I region where tensor B resides and where 
B| cns matters. This is what has been done in [la |. Un- 
der this condition Eq. © is a good estimate until the 
noise level is so low that the effects of neglecting the non- 
Gaussianity in iVj become present. 

As we pointed out earlier, this method obviously 
breaks down when one applies delensing to the search 
of B-mode at high I (such as cosmic strings or rotation 
induced -B-mode) . The biases are far from negligible even 
for a modest noise level. One could try to calculate all the 
bias terms and correct for them. However, it is difficult 
to analytically track all the bias terms in the quadratic 
estimator. The situation is even more complicated for the 
iterative estimator because of all the N-point functions 
introduced by iterating. Even if an expression is found 
for all the correction terms, their numerical evaluations 
can also be non-trivial. In the next section, we propose 
a revised delensing process that works at high I. 



V. A BIAS-FREE DELENSING PROCEDURE 

We choose not to remove the delensing biases by ex- 
plicitly calculating all the terms. Instead, we adopt a 
similar method used in [l6j . in which the measured B- 
map is high-pass filtered before the deflection field recon- 
struction. In [16| the i-cut was chosen to be at I ~ 150, 
leaving most of the lensing information intact. The cutoff 
ensures the instrumental noise at low I and the degree- 
scale primordial B-mode are not artificially "delensed" 
due to the biases discussed in the previous section (Fig 
Realizing that, we can revise this method for high I 
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FIG. 4. The simulation result with cutoff at / = 144 
(black vertical solid line) for r — 0.05 and r — 0. At 
low I, the residual power recovers the theoretical prediction 
(sum of instrumental noise and reconstruction noise, bias- 
free), while at high I there are contributions from terms 
such as <£° bs .E obs <^) c and (B noiBC * (l)B lcnB (l)) for r = 0, or 
(B pri *(l)B lcnB (l)) for r = 0.05. 



region. We proposed a modified estimator B dcl (l) instead 
ofEq. Q: 



B 



v (h) 

del 



(27T) 



B acl (l) = B ohs (l) 



2 /(Ma) h 



C EE 



N EE cf 



N 



(0 



(9) 



First, the Fourier space is divided into several non- 
overlapping annuli. The annuli are labeled by index 
N = 1,2.... The area of the annulus determines how 
many modes we lose in lens reconstruction. Instead 
of high-pass filtering, we implement a "notch filter" to 
the Fourier B-map for each region N. In other words, 
the modified estimator of lensing potential in region N, 
</>tv(0j would be the ordinary quadratic or iterative es- 
timators except that the modes located inside the an- 
nulus N are set to zero. As usual, the estimator for 
lensed B-mode, B l £ ns (l), will be calculated using 4>n(1), 
except that the result of B l £ ns (l) outside the annulus N 
will be ignored because they have been contaminated by 
delensing bias. Repeat the process for a different annu- 
lus, eventually we obtain uncontaminated B]^ ns (Z) cor- 
responding to each annulus. The last step is straight- 
forward. We compute the delensed B-mode B del (l) by 
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subtracting the combination of these uncontaminated es- 
timated lensed B-modes from the observed £?-mode. The 
estimated power spectrum from B dcl (l) is now free of the 
delensing biases. In Fig|5j we show that this procedure 
avoids the delensing biases and the delensing residuals 
now agree with theoretical expectations very well. 

Ignoring the B-map in the annulus where delensing 
takes place will obviously degrade the signal-to-noise. 
The reduction in signal-to-noise can be minimized by di- 
viding the B-map into smaller annuli, if one is willing 
to pay the price of proportionally increased computing 
time. This is still considerably simpler compared to a 
more analytic treatment of {B lcns *(l)B lcns (l)), which will 
have to be computed numerically anyway. Another sig- 
nificant advantage is, it is straightforward to implement 
this "moving notch filter" approach to iterative estima- 
tors to have the lowest possible delensing residuals (Fig. 

It is known that the cut sky effects make the neigh- 
boring modes in Fourier space correlated. In this case 
the removal of E, B modes in certain annulus described 
above would not guarantee perfectly avoidance of delens- 
ing biases. We might improve the behavior by extending 
these annuli to overlap (each side to AZ ~ 360/0) but pre- 
serve the modes B l ° ns {l) in the original annulus. We see 
this as yet another issue related to polarization delens- 
ing with map boundaries (in addition to other unresolved 
issues) and defer a full treatment to future work. 

With the "moving notch filter" approach, the residual 
is now free of the correlation terms that are plaguing the 
delensing process. However, to estimate the primordial 
B power, the residual power still needs to be subtracted 
from the delensed B-power estimator. If primordial B- 
mode exists, it will contribute to lens reconstruction error 
N? . In other words, the residual power is understand- 
ably a function of the primordial B power. This is clearly 
seen in Fig. [SJ This issue can be easily resolved with a 
standard iterative process. We start by assuming there 
is no primordial B and subtract the appropriate residual 
from the delensed B power estimator. If the estimated 
primordial B agrees with zero all is good. If not, we 
iterate the formula below until it converges: 

Cf p " = {B ldcl * (l)B ldc \l)) - Cf rra (Cf p "), (10) 

where B' del (l) is the same as Eq. ©, except that we alter 
the weighting in lens reconstruction by including a pri- 
mordial B-mode Cf pri . The second term Cf 168 (Cf pri ) 
is Eq. © except we replace Cf lcna by Cf '™ B + Cf prl 
in Eq. ([3]) when evaluating iV, . cf 1 " 11 ' remains the 
same because primordial -E-mode is much smaller than 
lensed i?-mode. In actual data analysis we need to have a 
precise prediction for "residual" and understand the con- 
vergence properties. In Fig. [5] we demonstrated percent 
level agreements between theory and simulations. 




FIG. 5. This figure is the main result of this paper. It shows 
the comparison of theoretical predictions and simulated de- 
lensing residuals using a joint EE and EB iterative esti- 
mator. With the "moving notch filtering" described in the 
text, the simulated delensing residuals are in excellent agree- 
ment with the theoretical predictions calculated from Eq. ((6]) 
without the biases described in section IV. For demonstra- 
tive purposes, we assume a string induced primordial B-mode 
G/i — 1CP 7 ). When the primordial B-mode is non-zero, we 

replace cf 1 "" 1 " in Eq. (|6} with cf lalls + C\ prl to predict the 
theoretical residuals. The annuli are chosen to be I — — 180 
and I — 180 — 1044 with spacing equals 108. Each simulated 
point contains 25 CMB realizations. 



VI. APPLICATION OF DELENSING FOR 
TENSOR AND COSMIC STRING GENERATED 
B-MODE 

For a given experimental survey speed it is important 
to optimize the survey width to obtain the best possible 
sensitivity on the target signal. We provide these "ex- 
perimentalists' guides" for both the tensor induced and 
cosmic string induced B-mode in Fig. [S] Plots simi- 
lar to these exist in various forms for tensor-induced B- 
polarization, quantified by tensor-to-scalar ratio r 9] [26] 
[27| But to our knowledge, iterative estimators have 
never been considered in these plots, so we produce them 
here for comparison. Astronomical foregrounds are also 
neglected as they can be removed with multi-frequency 
observations. However, effects of cut sky have been ig- 
nored so these plots still should not be regarded as final. 

In previous section, we showed that it is possible to 
delens at high I region. The delensing algorithm for- 
mally only applies to Gaussian primordial signals. The 
rotation-generated S-mode is Gaussian. Strictly speak- 
ing, both the string generated B-mode and the patchy 
reionization B-mode are non-Gaussian. As has been 
pointed out by 0], it is possible to simultaneously re- 
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construct the <j> field and the optical depth r field. That 
would be compatible with the moving notch filtering ap- 
proach proposed in the previous section. There is not 
enough information on the non-Gaussianity of string in- 
duced B-mode. For demonstrative purpose only, we ap- 
ply that method to the search of string-induced B-mode 
by assuming Gaussianity. Although cosmic strings do 
not source the majority of initial perturbations, they are 
still generic predictions of particle physics and cosmology 
so should be searched for in the data. Some theoretical 
models such as KKLMMT [29[ can induce inflation along 
with cosmic strings. Even the contributions from strings 
are severely limited by CMB temperature power spec- 
trum, it has been shown in simulation that string gener- 
ated B-mode as large as the lensing induced B-mode is 
still allowed by data [30| (see FigO}. 

For the stringy B-mode we assume the cosmic strings 
are smooth (a = 1) with segment RMS velocity 
v r = 0.65 and initial correlation length divided by ini- 
tial conformal time as 0.13. 

When we fix the parameters other than G/x, the power 
spectrum of B-mode will be proportional to (G/x) 2 so 
that we choose (G//) 2 to be our parameter for forecast. 
We use Fisher information matrix to be our estimate of 
ability to constrain (G^i) 2 : 



dC, 



BB , string 
I 



d{Gnf 



-)V 



BB\ 



(11) 



where C(G/j) 2 is Fisher limit for detecting quantity (G/i) 2 
(also the smallest detectable (G/x) 2 ), af B is the stan- 
dard error with which each multipole moment C BB can 
be determined. Summing over I > 0/180 makes sure 
no information is coming from modes with wavelengths 
larger than the survey size (here we neglect E-B mixing 
which will decrease the number of useful modes). We 
assume the delensed field will share the same covariance 
properties with a Gaussian field. Since we intend to do 
null test no string B-mode is included in of B . Therefore, 



BB 



V/ sky (2Z + l) ( ' ^ } ' 



(12) 

where / s k y is the fraction of sky we survey, lo 1 = 
4:Tt(s/Tcmb) 2 Aobs 25|. Quantity s is the detector sensi- 
tivity in units /xKsec 1 / 2 , a is the FWHM of the beam 
and tabs is the total observation time. Here we set 
s = 13/xKsec 1 / 2 and t b s = 1 year (so that when 
/ sky ~ 0.0024 or 6 = 10°, A p ~ 1.41/xKarcmin). 
The forecast results are presented in FigjBJ 
Tensor Induced B-mode When lensing is ignored, 
the r vs. width curve recovers early predictions by [26| . 
in which a survey width as small as ~ 5°seems to be 
favored. When lensing is treated as noise and no de- 
lensing is performed, the optimal survey width increases. 
The minimum near or above 100 degrees is caused by the 
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FIG. 6. The simple forecast for detecting tensor B-mode and 
cosmic strings. The horizontal axis is the sky coverage 6 ~ 

1 /2 

203/ sk ' y . The vertical axis is smallest detectable tensor-to- 
scalar ratio r (top) and strings tension G/i (bottom). We set 
s = 13pKsec 1/2 and s = l/iKsec 1/2 , 4 aremfnute FWHM, 
and observation time t Q bs = 1 year. 



tensor reionization peak near I ~ 6 (Fig. [T]). This ap- 
parent minimum drives some suborbital experiments to 
cover as much sky as possible at the risk of seeing the 
Galactic foreground emission. With quadratic, and espe- 
cially with iterative delensing, deep surveys become more 
and more favorable as the instrumental noise decreases, 
although it never quite recovers the "no lensing" curve. 
Non-Gaussianity at low noise level and E-B leakage from 
cut-sky can further prevent the r limit from reaching the 
predictions for iterative delensing. In addition to these 
factors, the choice of the survey width will be a prac- 
tical balance between the required knowledge on noise 
properties and foregrounds (especially for wide surveys) 
and polarization/beam systematics (especially for deep 
surveys). 
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String induced B-mode The peak at I ~ 600 in the 
stringy £?-mode spectrum (Fig. [I} drives the optimal sky 
coverage to 9 < 20° even without any delensing. With de- 
lensing, the optimal coverage reduces further to 9 ~ 10°. 
In this case, the lens cleaning can improve the detectable 
Gfi ~ 2.8 x 10" 8 to Gfi ~ 1.6 x 10~ 8 with 9 ~ 10° 
(s = 13/^Ksec 1 / 2 ), G/i~8x lO" 9 to Gfi ~ 4 x 10~ 9 with 
9 ~ 100° (s = 1/^Ksec 1 / 2 ) . 

VII. CONCLUSION 

In contrary to previous assumptions in the literature, 
we show in this paper that many high order correlation 
terms introduced during even the simplest delensing pro- 
cess are far from negligible. This can be avoided by 
spatial filtering the observed B-map prior to deflection 
field reconstruction, a practice that has been adopted for 
tensor- induced B-mode in 16]. We proposed a more 
comprehensive fix to the problem. With the new filtcr- 
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